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ENHANCED TIME OF ARRIVAL ESTIMATION USING 
REDUCED COMPLEXITY OPTIMAL PROCESSING 



BACKGROUND OF THE INVENTION 
5 The present invention generally relates to wireless signal reception, and particularly 

relates to higher resolution techniques for estimating a signal's time-of-arrival. 

Position determination finds broad utility across a range of applications, including 
emergency 911 (E-911) locating services, navigation systems, location-based information or 
marketing services, and others. Determining the position of a wireless receiver typically 
1A involves some type of triangulation approach involving time-of-arrival (TOA), direction-of-arrival 
,g (DOA), or other like technique, wherein the receiver determines appropriate parameters for 
signals received from three or more transmitters. 

For example, a Global Positioning System (GPS) receiver determines the TOA for 
signals received from four or more GPS satellites to uniquely determine the receiver's latitude, 
f§ longitude, and elevation. Wireless access terminals, such as the mobile radiotelephones used 
q in modern wireless communication networks, may also derive their position through triangulation 
p based on determining the TOA of signals received from three or more geographically separated 
radio base stations. Conversely, determining the TOA of the receiver's signal at each of the 
three radio base stations yields the same result. In any case, the accuracy of position 
20 estimation depends on the accuracy of TOA estimations. 

Accurate TOA estimation entails a number of challenges. Some approaches to TOA 
estimation, while providing accurate, high-resolution estimates, are computationally too complex 
to yield accurate position estimates within the requisite times. For example, emergency locating 
requirements applied by the Federal Communication Commission to wireless network operators 
25 for 91 1 calls requires relatively accurate position determination (50-100 meters or less) within 



m 
in 



1 



Ericss^^^l 



if. No. P12633-US1 RCUR 
C&BRef. No. 4015-987 



five seconds or less. In these circumstances, using highly accurate but potentially very slow 
TOA estimation algorithms is unacceptable. 

Other challenges arise from the real-world vagaries of the typical received signal. What 
originates as a single transmit signal oftentimes arrives as a composite received signal 
5 comprising a plurality of multipath rays. Each multipath ray travels a different propagation path, 
and thus has different signal characteristics in terms of phase shift and attenuation, and, 
significantly, each ray has a different TOA with respect to the wireless receiver. 

Thus, estimating TOA for the received signal entails processing the multipath signal 
such that at least the significant multipath components are resolved from the composite 
1© received signal. In some cases, the time separation between rays will be minimal, further 
^ complicating the task of resolving the multipath rays with enough resolution to satisfy the overall 
^ TOA estimation accuracy requirements of the positioning task at hand. 



l, BRIEF SUMMARY OF THE INVENTION 

□ 

ijK Received signal TOA estimation in accordance with at least some embodiments of the 



present invention combines sub-optimal and optimal received signal processing to approach the 
accuracy obtainable only with a pure optimal approach, while simultaneously reducing the 
aggregate computational complexity. Pre-processing the received signal with a sub-optimal 
algorithm reduces the search space within which the computationally more intensive optimal 
20 algorithm, such as Maximum Likelihood Estimation (MLE), is applied. The extent to which the 
search space is reduced or decreased in dimensionality depends on the particular sub-optimal 
approach taken. In general, sub-optimal processing restricts at least the time window or 
windows within which the optimal algorithm operates. 

Thus, sub-optimal processing may be used to simply identify one or more time intervals 
25 over which the higher resolution optimal algorithm is applied. Depending upon the time 
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separation of the multipath rays, sub-optimal processing may be used to identify the probable 
number of multipath rays in the received signal, and further to generate rough TOA estimates 
for those rays. Optimal processing may then operate at higher resolution within time windows 
defined around the estimated ray positions. MLE processing operates within these time 

5 windows to identify the most likely combination of ray arrival time by considering all 

permutations of ray position, amplitude, and phase for the one or more multipath rays involved. 
The position of the earliest ray corresponding to the Maximum Likelihood permutation may be 
taken as the TOA of the received signal. 

Sub-optimal processing of the received signal may entail using any sub-optimal 
1§3 algorithm, with correlation, Multiple Signal Identification and Classification (MUSIC), and Signal 

h *0 Eigen Vector (SEV) approaches representing exemplary algorithms. In some cases, a sub- 

y optimal algorithm is used iteratively to further refine TOA information in advance of optimal 

III 

;^ algorithm processing. For example, correlation processing may be performed iteratively, with 

^ successive cancellation of correlation peaks in each iteration cycle. This approach iteratively 

k-J 

# breaks the received signal down into its ray components, and may, in some instances, provide 

W 

2 sufficient accuracy for TOA estimation absent any subsequent optimal processing. Of course, 

^ optimal processing may still be applied within time windows defined around one or more of the 
ray components. 



20 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram of a wireless device receiving a transmitted signal as a composite of 

several multipath signals. 

Fig. 2A is a diagram of multipath components in a received signal with relatively wide 

time separation. 
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Fig. 2B is a diagram of multipath components in a received signal with relatively narrow 
time separation. 

Fig. 3 is a flow diagram of an exemplary approach to time of arrival estimation using a 
combination of sub-optimal and optimal algorithms. 
5 Figs. 4A-4C are diagrams illustrating the successively better resolution of correlation 

peaks in the received signal using iterative sup-optimal processing. 

Fig. 5 is a flow diagram of iterative sub-optimal processing. 



1© DETAILED DESCRIPTION OF THE INVENTION 

<£ 

^ The present invention broadly applies to received signal TOA estimation. It is not limited 

to use in mobile terminal positioning within wireless communication networks, although the 
invention may be used to great advantage in such environments. Broadly, TOA estimation in 
accordance with the present invention applies to a wide range of signal types and systems. 



In general, the present invention uses a sub-optimal algorithm to reduce or otherwise 



^ restrict the search space of an optimal algorithm used to precisely estimate the arrival time of a 

y 

received signal. Initial processing with the sub-optimal algorithm reduces the optimal 
algorithm's search space by generating crude estimates for one or more multipath components 
of the received signal. In at least one embodiment, the sub-optimal algorithm permits 
20 identification of a number of probable multipath components of the received signal, such that the 
optimal algorithm may be restricted to time intervals around one or more of these probable 
multipath components rather than across a more expansive time window. This method allows 
the optimal algorithm to operate at fine time resolutions without incurring impractical 
computational complexity. 
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Turning now to the drawings, Fig. 1 is a diagram of a typical radio signal transmission 
between a transmitter 20 and a receiver 30. Because of multiple propagation paths, denoted as 
Li through l 3 , the transmitted signal from transmitter 20 is received at receiver 30 as a composite 

of multipath signals or rays, with each ray traveling through one of the propagation paths l. This 
5 composite received signal may be denoted as r'ft). Typically, the receiver 30 filters, digitizes, 
and down-converts the initial received signal r'(t) to create the received signal rft), which is a 
composite of the multipath signals. That is, rft) includes multiple versions of the transmitted 
signal s(t), with each version having its own characteristic channel distortions affecting signal 
phase, amplitude, and time of arrival. In general, these multipath components complicate 



II] particular application in which it is used. For example, for use in wireless communication 

m 

U systems, the receiver 30 may be configured in accordance with TIA/EIA/IS-136 standards, 

g Tl A/El A/I S-95B standards, TIA/EIA/IS-2000 standards, or any of the other common air interface 

fS 

1il standards, such as the European GSM and Wideband CDMA (W-CDMA) standards. Generally, 



0 the receciver 30 will comprise both optimal and sub-optimal processors, which may themselves 
comprise all or part of one or more digital signal processors (DSP) within the receiver 30. 

Fig. 2A illustrates an exemplary signal rft) comprising four apparent multipath signal 
components corresponding to arrival times 7"? through T 4 relative to a reference time T 0 . Here, 

20 the multipath signals appear to be widely spaced relative to the symbol timing used to carry 
information in the received signal, and sub-optimal processing enables identification of the four 
probable multipath components. Even so, each multipath component may in reality comprise 
two or more closely spaced multipath signals, which can be accurately resolved by applying the 
optimal algorithm to time windows defined around each one of the four probable multipath 

25 signals identified by sub-optimal processing. 




processing the received signal rft) to estimate its time of arrival. 



In the context of the present invention, the receiver 30 is configured based on the 
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References Wi through W 4 denote the time intervals or windows defined around times Tj 
through T4. As will be explained in more detail later, optimal processing is restricted to time 
windows W } through W 4 in the interest of reducing the optimal algorithm's search space. Using 
the optimal algorithm, the precise arrival time of the received signal r(t) within the time windows 
Wi through W 4 may be more accurately fixed relative to the reference time T 0 . Further, if one or 
more of the four apparent multipath signals is itself comprised of two or more closely spaced 
multipath signals not resolvable using sub-optimal processing, optimal processing will resolve 
the additional multipath signals within its time resolution limits, thereby allowing precise TOA 
estimation for one or more multipath signals within each of the time windows W } through W 4 . 

Fig. 2B illustrates another exemplary signal r(t) f but rather than having relatively widely 
spaced multipath signal components, its multipath components are essentially overlapping with 
respect to the symbol timing of r(t). Indeed, this illustration may depict the scenario described 
above where one of the four apparent multipath signals shown in Fig. 2A itself comprises two or 
more closely spaced multipath signals. Where the degree of overlap is great, sub-optimal 
processing may not identify more than one general time window over which the optimal 
algorithm is used. Here, W } denotes the time window spanning the closely grouped multipath 
signal components of r(t). The time window W l may be formed by applying a threshold function 
to the composite received signal rft), illustrated by the dashed-line pulse superimposed over the 
group of multipath signal components. Other techniques, such as least squares processing 
may also be used to define the span of the search windows W. 

Thus, with at least some embodiments of the present invention, TOA estimation involves 
at least a two-step process. First, a sub-optim al algorithm generates crude est imates of arrival 
times of multipath components in the received signal r(t). Second, an o ptimal algo rithm r efine s 
these crude estimates, such as by using Maximum Likelihood (ML) techniques, although in 
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some cases a near optimal algorithm such as SEV may be used. Search spaces for the 
multipath signals in the second step are restricted to within the vicinity of the crude estimates. 

While essentially any sub-optimal approach may be used to initially process the received 
signal, exemplary sub-optimal algorithms include correlation, MUSIC, and SEV. While these 
5 algorithms are subject to variation, exemplary details regarding their use with the present 

invention are provided below. However, it may be helpful to begin with a brief description of an 
exemplary mathematical model of the received signal. 

Let Q copies of the transmitted signal be received by the receiver 30, where the q th copy 
of the received signal is given as, 

If r<*>«= Za]Sf^-T m )+P^(0, qe{l 9 2,... 9 Q} (1) 

%Q m=l 

where s(t) is the transmitted signal, x m is the time delay of the m th path, affi is the 

m 

M, attenuation of the m th path due to scatter characteristics and propagation through the medium, 

O and p^(/)is additive white Gaussian noise (AWGN). In the above model, the succession of 

CO 

□ received signals r Kq) (i) is considered as resulting from repeated independent experiments on 

1§ the medium. The time delays of the paths x m may be considered to remain unchanged for each 

of the Q received signals. Also, affi and a^f * may be considered to be uncorrelated when 
q * p . The problem to be solved is to obtain the estimates of the path delaysx w . With 
accurate estimates of the path delays, the TOA of the composite received signal may be 
accurately estimated as, for example, the TOA of its earliest multipath signal component. 
20 Typically the received sequences are sampled so that the model considered in (1) 

becomes, 

r<*>(n)= fa^s(n-x m ) + ^\nl q e {1,2,..., 0},/? e {I..N}. (2) 
m=l 
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That is, each received sequence is N samples long. In matrix form, this may be 
represented as, 

rM=S-*M+bte\ (3) 
where and 

„(<?) 

are column vectors of length N t a^ is a column vector of length M, 

and, 



S = 



S(l-Tl) 5(1-T 2 ) • • • S{\~X M ) 

s(2-x { ) 5(2 -x 2 ) • • • s(2-x M ) 



s(N-xi) s(N-x 2 ) • - • s(l-x M ) 



(4) 



The elements of are the channel attenuation valuesa^ where m e {l.M} , and the 

elements of are the noise values p^(«) , where n e {I..N} . The columns of the matrix S are 
the signal vectors corresponding to each delayed multipath component of the received signal 

r(t). 

Correlation processing represents an exemplary approach to sub-optimal processing of 
the received signal r(t). With correlation processing, the received signal r(t) is correlated with 
the transmitted signal corresponding to each one of a number of delay hypotheses, and the 
delay hypotheses that yield maxima in the energy of the correlator output are chosen as the 
path times of arrival. Thus, the correlator output may be described as, 



c(x) = 5(x) // r, 



(5) 



where s(x) - [s(l - x) s(2 - x) . . . s(N- x)y, [•]" indicates the Hermitian matrix, and where 
[•] T indicates the transpose matrix. Correlation techniques and apparatus that might be included 
within the receiver 30 for correlating the received signal r(t) are well understood by those 
ordinarily skilled in the art. 
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With regard to wireless communication applications, the received signal r(t) may be in 
essentially any form, as in Time Division Multiple Access (TDMA) systems based on, for 
example the TIA/EIA/IS-136 standard, or as in Code Division Multiple Access (CDMA) systems 
based on, for example, the Tl A/El A/I S-2000 standard. With these air interface definitions, and 
5 with others, such as the European Wideband CDMA (W-CDMA) standard, the received signal 
r(t), or a signal associated with the received signal, contains known information, which may be 
used by the receiver 30 for correlation operations. Such known information may comprise any 
known sequence of data or symbols, such as a training sequence or pilot symbols, or may 
comprise a known spreading code. In any case, the receiver 30 may base its correlation 

1p operations on this known information, and, as described earlier, may accumulate multiple copies 

q 

^ (e.g., Q copies) of this known information by periodically receiving known information from the 
UJ 



m 



ill 



transmitter 20. 

A modified form of the above correlation operation may be used when the correlation 
matrix K N of the impairment process $(n) is known to not be an identity matrix, i.e., the 



u 

impairment process is correlated. This form of the correlator output is given by, 



c(z) = (K^)) H r, (6) 

where K^ 1 is the inverse of the noise correlation matrix K N . The noise correlation matrix 

K N in general indicates the correlation of the noise between various received samples. If the 
noise is uncorrelated, the matrix has non-zero values on its main diagonal only, as noted above. 
20 The correlation matrix of the received signal sequence (samples of r(t)) is given by, 

K = E(r-r H ) = SPS H +K N , (7) 
where P is the correlation matrix of the received multipath components of r(t) whose 

elements are given by Py = £(<x; -a*) . The matrix K is estimated asK by averaging the Q 

received copies of the received signal. Thus, the estimated correlation matrix is given as, 

9 
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(8) 



As an alternative to the exemplary correlation process details, sub-optimal processing of 
the received signal may employ the MUSIC algorithm, which entails finding the Eigen values 

and Eigen vectors of the matrix K as follows, 



where the columns of V contain the Eigen vectors oik and D is a diagonal matrix given 
byD = diag(\ h ...,X N ) , where Xi are the corresponding Eigen values. The MUSIC algorithm thus 
estimates the delays of the M paths for of the multipath components of the received signal r(t) 
as the M maxima of the function, 



where V N is the matrix containing the TV- M columns of Fthat correspond to the noise 
Eigen vectors of F, and s(z) is as defined earlier in (5). The noise Eigen vectors of Fare the 
columns of Fthat are orthogonal to the signal vectors that make up the received signal r(t). The 
delays of the multipath signal components in r(t) are estimated as the delays x, that correspond 
to the peaks of the function given above in (10). 

The MUSIC algorithm may thus be summarized by the following steps: 

(a) receive Q copies of the sampled received signal r(t) given by (2); 

(b) form an estimate £ of the correlation matrix K as given by (8); 

(c) find the Eigen values from the matrix D and the Eigen vectors Fof the matrix 

K according to (9); 

(d) estimate the number of multipath components M from the Eigen values in the 
diagonal of the matrix D, where the high Eigen values correspond to the signal 
components; 



KV - KjsjVD , 



(9) 




(10) 
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(e) collect a set of noise Eigen vectors in a matrix V Ni where the noise Eigen vectors 
are the Eigen vectors in Fthat correspond to the low Eigen values in the diagonal 
elements of D\ 

(f) for each hypothesis of a delay value, x e(T min ..T max ), compute the value of the 
function F M usic(^) as given by (10); and 

(g) find the delay hypotheses x m for which the function Fmjsi&x) achieves its M local 
maximum values, where the ^corresponding delays(T l5 T 2 ,. ..,T^)are the 

estimates of the times of arrival for the multipath components of r(t). 
As the MUSIC algorithm involves finding the M maxima of a one-dimensional function 
1 <j instead of the maximum of an M-dimensional function, the MUSIC algorithm has significantly 

l^j lower complexity than that of the optimal, Maximum Likelihood approach (ML). 

Wj 

jji The SEV approach involves the maximization of the M-dimensional function, 

L F SEV (x h ...,x M ) = TrOloKxVsWvi 1 ), (1 1 ) 

g where n 0 = (/ - U s ) . with U s = S^~ l S H K^ and ¥ = S H K^S , I is an identify matrix, V s 

1|5 is a matrix whose columns contain the columns from the matrix Kthat correspond to the M 

signal components of r(t)X^, the signal Eigen vectors, and Wis a weighting matrix. It may be 

shown that the optimal weighting matrix is given by, 

W = [D m -I] 2 Dm\ (12) 
with D M ^diag{X\,...,X M ) , where "diag" indicates a diagonal matrix, and where A,, are the 

20 Eigen values corresponding to the signal Eigen vectors in consideration. 

The SEV approach may be summarized in the following steps: 

(a) receive Q copies of the sampled received signal r(t) as given by (2); 

(b) form an estimate k of the correlation matrix K as given by (8); 
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(c) find the Eigen values from the matrix D and the Eigen vectors Fof the matrix 
K according to (9); 

(d) collect a set of signal Eigen vectors in a matrix V s , where the signal Eigen 
vectors are the Eigen vectors in V that correspond to the M high Eigen-values in 
the diagonal elements of D\ 

(e) compute the weighting matrix Was given by (12); 

(f) for each hypothesis of the delays i M of the M paths, compute the values 

of FsEvb\>*2>->*M)as 9 iven by (11); and 

(g) find the hypothesis for which the function F SE v achieves its maximum value, 
where the corresponding delays (ti,t 2 v ) x m) are the estimates of the times of 
arrival of the M multipath components of r(t). 

In the SEV approach, an estimate of the number of multipath components may be made 
from the Eigen values. However, where M is not known, it must be estimated to set the 
dimension of the F SE v function, which must be maximized over the M estimated paths. The 
complexity of the SEV algorithm may approach that of the Maximum Likelihood (ML) optimal 
approach, which is detailed below. 

With the multistage TOA estimation techniques used in at least some embodiments of 
the present invention, sub-optimal processing, as exemplified above, reduces the dimensionality 
of the search space within which optimal (e.g., ML) processing is applied. This yields the TOA 
estimation accuracy close to that obtainable only with optimal algorithms, yet reduces the 
aggregate complexity of the overall TOA estimation process. 

A more detailed look at optimal processing based on the ML approach begins with a 

definition of the ML function. After defining X ¥ = S H K~ { S , the ML function is given as, 



F ML ( T 1 > T 2 v-,x M ) = Tr(Tl s kK N ] ) , 



(13) 
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where Tr(.) is the trace operator and n s =s^- ] s H K~ l , which is an Nx N projection 

matrix onto the subspace spanned by the signal vectors in the columns of the matrix S. 

This multistage sub-optimal/optimal approach may be summarized in the following steps: 

(a) process the received signal using a relatively low-complexity algorithm (as 

5 compared to ML processing complexity) to obtain initial estimates of arrival times for one or 
more probable multipath components of the received signal, with exemplary techniques 
including correlation and MUSIC processing as outlined above; and 

(b) apply optimal processing such as ML based TOA estimation within the reduced 
search space or spaces corresponding to the initial delay estimates identified during the initial 

a 

1CB processing. 



Viewed another way, optimal processing compares estimates of the received signal 



In against what was actually received, where the estimate of the received signal is based on 

111 

M= combinations of the variables affecting the received signal. The combination of variables 

O resulting in the estimate with the lowest error is taken as the most correct. That is, the estimate 

S3 

1^ corresponding to that set or combination of variables represents the "maximum likelihood" 



among all of the estimates. 

Because ML processing considers all combinations of variables, its complexity increases 
rapidly with the number of variables involved. Assume that the number M of multipath signal 
components in a received signal is known. The ML processor must compute an estimate of the 
20 received signal for all possible combinations of time delay, phase, and magnitude for all M 

multipath signals. If the number of multipath signals is not known, the ML processor must also 
consider all possible values of M, or at least must consider all possible number of multipath 
signals within a defined upper limit for M. 

Obviously, any pre-processing that restricts this variable space can result in significant 
25 reductions in the total number of optimal processing calculations that must be performed. 
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Within the context of the present invention, sub-optimal processing can, for example, identify the 
number of likely delay components within the received signal, and generate rough or initial time 
delay estimates for these identified components. 

Hence, the search space for optimal processing and thus computational complexity may 
5 be reduced by restricting the time windows for the initial delay estimates x m to a window of time 

defined byi m ±5 m , where 5 m is much smaller than the spacing between the delay estimates f m . 
This approach is illustrated in Fig. 2A, where the time windows W1 through W4 are defined 
around the estimated multipath delays of the received signal r(t). Here, the parameter 5 m 
defines the size of the ML search window around the vicinity of each delay estimatef m . Each 

O 

1C0 delay estimate-^ may use a different value for the 5 m parameter, and the span of the search 



window may be further adjusted based on the magnitude of the output and the shape of the 
output in the vicinity of the local maximum for the estimated m th path of the received signal. 
Where two time delay estimates satisfy (cV^ )<t , where / is a defined time difference 

CO threshold value, optimal algorithm search spaces may be modified through hypotheses 

W 

1£g including a third multipath component with a delayf* such that three paths within a delay 

0 

r* window (f ; -8 z y) <i k < (ij +8,y) are tested during the optimal processing phase. That is, 

hypotheses within this range are tested using an ML or possibly a SEV approach to further 
refine the TOA estimates. Thus, when two multipath signal components occur very close to 
each other in time within the received signal, a third multipath component is considered in the 
20 vicinity of the two closely space multipath components. The quantity 5, 7 is a parameter that 
defines the time span of the search window that will be used to search for the presumed third 
path within the vicinity of the closely spaced i th and j th paths. 

Where an estimate of the slope of the function F MUSIC (i.e., equation (10)) is low in the 
vicinity of a peak at a time delay i m , this indicates a multipath component in the received signal 
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r(t). At such points, optimal algorithm search spaces may be modified by hypothesizing second 
multipath components within a delay given asf m ±5 m . This allows the optimal search algorithm 
to search for additional multipath components within a restricted search space around the 
crudely identified multipath components revealed during sub-optimal processing. 

Thus, with the time window restriction, the ML algorithm may be set to test hypotheses 
of the various combinations of multipath signal component delay, amplitude, and phase, across 
a set of defined points within the search window or windows. The resolution of the time 
positions within the search windows at which ML hypotheses are evaluated for error with 
respect to the actual received signal may be set depending upon the time resolution 
requirements placed on TOA estimation. Higher resolutions require subdividing the search 
windows for ML algorithm processing into finer subdivisions of time. For example, the ML 
algorithm may be set to test hypotheses within one or more search windows at time positions 
spaced apart at one-tenth the received signal's basic symbol timing for accurate TOA 
estimation. Where less accuracy is required, the time positions may be set more coarsely. 

Fig. 3 illustrates exemplary flow logic for combining sub-optimal and optimal algorithms 
in TOA estimation. Processing begins with reception of the received signal r'(t), which is 
converted to baseband to form the sampled received signal r(t) (step 300). The sampled 
received signal is the processed using, for example, the MUSIC algorithm or the correlation 
method to identify one or more multipath signal components of the received signal (step 302). 
This pre-processing of the received signal provides crude estimates of the actual multipath 
components in the received signal, and serves to reduce the search space used in optimal 
algorithm processing. Here, sub-optimal processing has produced M path estimates. 

Next, the path delays of the multipath signal components identified in sub-optimal 
processing are sorted in ascending order (step 304). That is, the earliest TOA is placed first in 
the sort order. A path index variable / is then set to the first path in the sort order (step 306). 
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Path pairs that are adjacent in time based on the sort order are then evaluated to determine 
their proximity to each other. This is in keeping with the earlier point regarding hypothesizing 
additional multipath components within the vicinity of two closely spaced multipath components 
identified during sub-optimal processing. Specifically, the time difference between the i th 
5 multipath component delay and the next multipath component delay in the sort order (i.e., the i th 
+ 1 or j th delay) is evaluated to determine whether it is less than the defined time difference 
threshold /(step 308). 

If the time difference between the i th multipath component delay and the next multipath 
component delay in the sort order exceeds the threshold t, the optimal algorithm search window 
l!3 Wi for the i th multipath component delay is defined as f ; -5/ \o% { +5 i9 whereSj is a defined time 
|| span that may be set or adjusted on a per-delay basis, depending on, for example, the slope of 
yi the output in the vicinity of the corresponding i th multipath component, and wheref j is the delay 
M relative to some reference time (e.g., T 0 ) of the i th component's estimated delay time (step 310). 

p The component index variable i is then incremented (step 312) and evaluated to determine if 

00 

1i5j there are additional multipath components as identified by sub-optimal processing (step 314). If 

y 

O not, optimal processing begins within the restricted search windows. If additional components 
remain, the window restriction process is applied to the delay time of the next component in the 
sort order (steps 308+). 

If, however, the time difference between the delay corresponding to the i th component 

20 and the i th + 1 component does not equal or exceed the threshold t (step 308), the j th component 
(j > 1) at which the corresponding time delay satisfies the equation f y -ij_ x > t is identified (step 

318). Essentially, this entails incrementing the index j from j = i + 1 until the time difference 
between the j th and Q - 1 ) th multipath component delays exceeds the defined threshold u Then, 
an optimal algorithm time window is defined that spans the ith through the jth time delays as 
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f -5 ; tof j +8 y , where t, is the delay time of the i th component and f j is the delay time of the } 



:th 



component (step 320). Note that this time interval includes one or more delay times for 
multipath components having delays greater than the i th delay but less than the j th delay. After 
defining this search window, the sort index i is incremented to one past the jth index determined 
in step 318 (step 322), and processing continues (step 31 4+). 

While combining sub-optimal and optimal processing yields the most accurate received 
signal TOA estimates, there may be some instances where the required accuracy in TOA 
estimation does not necessitate optimal algorithm processing. In these cases, sufficient 
accuracy may be obtained by iterative processing of the received signal using a sub-optimal 
approach. For example, the correlation method may be used iteratively in a successive 
cancellation approach. With successive cancellation, the multipath signal component, or 
components identified in each iteration are treated as noise and cancelled for the next iteration. 
With each cancellation, the actual multipath components of the received signal are revealed 
with successively better resolution. 

In more detail, at the beginning of each iteration a new noise correlation matrix K N is 

calculated a$K N = K N +S'PS' H , where the matrix S' has columns that are the signal 

components at the estimated time delays from the previous iteration. This effectively treats the 
detected signal vectors (i.e., the vectors corresponding to the detected multipath components of 
the received signal) as noise. The matrix P' is the correlation matrix of the channel for the 
corresponding received signal paths. If this matrix is not known, it may be estimated or chosen 
as the identify matrix. 

Then, for the next successive iteration, the sub-optimal algorithm, such as the MUSIC or 
correlation algorithm, is re-executed as previously described, with Aroused rather than K N . The 
output of each iteration may be used to determine if one or more additional iterations are 
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needed. For the correlation method, the energy at the output of the correlator is used to 
determine whether additional paths may be present. For the MUSIC algorithm, the sum of the 
Eigen values computed at each iteration is used to decide whether additional paths may be 
present. In either case, if the metric is above a certain threshold indicating that signal energy 
5 from more paths is present, another iteration is performed. If the metric falls below the specified 
threshold, further iterations are deemed unnecessary, and the multipath component detected 
that has the lowest path delay is chosen as the first path. The TOA estimate for this first path 
may be taken as the TOA estimate for the received signal. 

Figs. 4A-4C illustrate iterative correlation as applied to a received signal having six true 
1£J3 multipath components. The "x" markers along the top grid in the figures identify the true time 
*0 delays of the multipath delays relative to the zero point on the horizontal axis. Fig. 4A illustrates 
yJ the initial correlator output for the received signal, which indicates a single, rather broad 

in 

W correlation peak. This response is similar to that depicted in Fig. 2B, where the delay 

y> 

* separation between multipath components of the received signal were too small to resolve in a 
1§P single-pass sub-optimal approach. 

9 With the iterative approach however, a delay time corresponding to the center of the 

H* single correlation response peak is assumed, and the noise matrix K N is updated such that the 
presumed multipath signal component corresponding to this delay time is cancelled. Such 
cancellation reveals additional detail regarding the true nature of the multipath components in 
20 the received signal, as illustrated in Fig. 4B, which depicts the correlator output on the second 
iteration. With multipath component identified in the first iteration cancelled, this second 
iteration reveals three correlation peaks. The presumed multipath components corresponding 
to these three peaks are then cancelled from the received signal (e.g., treated as noise 
components) by adjusting the noise matrix K N used in correlation operations. 
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Fig. 4C illustrates the path delay estimates on the third iteration after two rounds of 
cancellation. Here, the approximations of the six true multipath components of the received 
signal are more clearly revealed, and are positioned at relative delay times (refer to the "x" 
markers) at or close to the actual delay times of the true multipath components. It may be that 
in some applications, these delay time approximations arrived at iteratively are accurate 
enough. In this case, the TOA estimate of the received signal may be, for example, taken as 
the arrival time of the earliest multipath component illustrated in Fig. 4C. If these TOA estimates 
are not sufficiently resolved for the purpose at hand, they may be further refined by optimally 
searching time windows about one or more of them. 

That is, optimal processing may be combined with iterative sub-optimal processing. 
Such a combination may have advantages where optimal processing will be used at high time 
resolutions, such as less than one-tenth of the symbol timing used in the received signal. At 
such high resolutions, optimal searching even in reduced search spaces is computationally 
intensive, and entering optimal processing with strictly reduced search windows is 
advantageous. 

Fig. 5 illustrates exemplary flow logic for iterative processing. Operations begin with 
converting the received signal r'ft) to baseband to form the received signal r(t) (step 330). 
Then, sub-optimal processing, such as by MUSIC or by correlation, is applied to the received 
signal to estimate the number of multipath components in r(t), and to generate TOA estimates 
corresponding to these estimated paths (step 332). 

The matrix S' is estimated using the signal vectors at the estimated delays (step 334), 

and a new impairment correlation matrix (noise matrix) is calculated as K N =K N +S'PS* H (step 
336). Then, the sub-optimal algorithm is applied in a next iteration using the new impairment 
correlation matrix a:^ to estimate multipath component path delays (step 338). If more paths 
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are detected in this current iteration (step 340), then processing repeats for the next iteration 
using the adjusted impairment correlation matrix K N (steps 334+). If no additional paths were 
detected on the current iteration (step 340), processing continues with aggregating path delays 
across all iterations to form a final set of estimated path delays (step 342). 

Thus, TOA estimation may be based on a sub-optimal/optimal approach, an iterative 
sub-optimal approach, and an iterative sub-optimal/optimal approach. The above discussion 
provides details regarding exemplary approaches for these various combinations of sub-optimal 
and optimal received signal processing. In general, the present invention reduces the 
complexity of highly accurate TOA estimation using optimal algorithm approaches by preceding 
optimal algorithm processing with sub-optimal algorithm processing. This two-step approach 
may also be used to reduce the complexity of the SEV-based approach to TOA estimation for 
the multipath components of the received signal. Further, in some instances, iterative sub- 
optimal processing may obviate the need for subsequent optimal processing where TOA 
estimate accuracy requirements are not too stringent. 

Obviously, the present invention is subject to much variation. Thus, the exemplary 
details above should not be construed as limiting the present invention. Rather the present 
invention is limited only by the scope of the claims below, and the reasonable equivalents 
thereof. 
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